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NONLINEAR SIMULATIONS OF THE CONVECTION-PULSATION COUPLING 

T. Gastine^ and B. Dintrans^ 

Abstract. In cold Cepheids close to the red edge of the classical instability strip, a strong coupling between 
the stellar pulsations and the surface convective motions occurs. This coupling is by now poorly described by 
1-D models of convection, the so-called "time-dependent convection models" (TDC). The intrinsic weakness 
of such models comes from the large number of unconstrained free parameters entering in the description 
of turbulent convection. A way to overcome these limits is to compute two-dimensional direct simulations 
(DNS), in which all the nonlinearities are correctly solved. Two-dimensional DNS of the convection-pulsation 
coupling are presented here. In an appropriate parameter regime, convective motions can actually quench 
the radial pulsations of the star, as suspected in Cepheids close to the red edge of the instability strip. These 
nonlinear simulations can also be used to determine the limits and the relevance of the TDC models. 
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1 Introduction 

The cold Cepheids located close to the red edge of the classica l instability strip have a lar^e surface convective 



zone that affects their pulsation properties (e.g. the reviews of lGautschv fc SaiolllQQd : lBuchlerll2QQ9[ ). The first 



calculations, that assumed frozen-in con vection, predicted a cooler red edge than the observed one. Indeed, as 
already stated bv lBaker fc KippenhahnI ([l965), a non-adiabatic treatment of the convection-pulsation coupling 
is mandatory to predict the red edge location with a better accuracy. 

Several time- dependent convection (TD C) models were therefore developed to address this coupling (e.g. 
IStellingwerdfTosit iKuhfuiBlfToset IXiondll989h and succeeded in reproducing the co rrect location of the red edge 



despite their disagreements w ith the physical origin of the mode stabilisation (e.g.l Bono et al. I I1999I : Yecko et al 
19981 : iGrigahcene et al.ll20Q5[ ). 



However, all these formulations involve many free and degenerate parameters (e.g. the seven dimensionless 
a coefficients used by lYecko et~al] Il998[ ) that are either fitted to the observations or hardly constrained by 



theoretical values. Nevertheless, another way to tackle this problem is to compute 2-D and 3-D direct numerical 
simulations (DNS) that correctly take into account the nonhnearities involved in this coupling. Results of such 
pioneering 2-D nonlinear simulations of the convection-pulsation coupling are presented in the following. 

2 The convection-pulsation coupling 

Our system corresponds to a local zoom around an ionisation region responsible for the driving of the acoustic 
modes excited by the /^-mechanism. It is composed of a 2-D cartesian layer filled with a monatomic and perfect 
gas. The opacity bump associated with this ionis ation zone is modelled by a temperature-dependent radiative 
conductivity profile K(T) (for further details, see lGastine fc Dintrans|[2QQ8al ). In addition to the a^- mechanism, 
the conductivity profile is deep enough to locally get a super adiabatic temperature gradient, meaning that 
convective motions develop there, according to Schwarzschild's criterion. A strong coupling between convection 
and the acoustic oscillations therefore develops. 

The hydrodynamical equations are then advanced in time with the high-order, finite-difference Pencil Cod^^ 
which is fully explicit except for the radiative diffusion term that is solved implicitly thanks to a parallel alternate 



^ Max-Planck-Institut fiir Sonnensystemforschung, Max-Planck-Strasse 2, 37191 Katlenburg-Lindau, Germany 
2 I RAP, CN RS-Universite de Toulouse, 14 avenue Edouard Belin, F-31400 Toulouse Fran ce 
* See |http : //w ww . nordita . org/sof tware/pencil- code/ and IBrandenburg DobleJ (|2002l ). 
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direction implicit (ADI) solver of our own ( Gastine fc Dintraiisll2QQ8a ). The simulation box spans about 10% 
of the star radius around the ionisation region. In order to ensure that both the thermal relaxation and 
nonlinear saturation of the A>:-mechanism are achieved, the simulation s are computed over more than 4000 days 
(corresponding roughly to 1500 periods of oscillation, see lGastine fc Dintranslboilah . 

Figure [1] displays a snapshot of the vorticity field for two simulations discussed in'Gastine fc Dintrans' ('2011a'), 
namely G8 (upper panel) and G8H8 (lower panel). This vorticity field highlights the convective motions that 
are approximately localised in the middle of the layer, where the radiative conductivity is minimum. Differences 
in the typical length-scale of convection are noticeable between these two DNS: convective eddies are smaller 
scale in the G8 simulation than in the G8H8 one. Accordingly, the overshooting of convective elements into the 
lower stably stratified layer is also more pronounced in the latter simulation. 

Beyond these qualitative differences, a good way to compare these simulations is to study the temporal 
evolution of average quantities, such as the vertical mass fiux pUz- Indeed, as we are considering simulations with 
both convective motions and oscillations of acoustic modes, it is relevant to use a simple diagnostic that roughly 
separates their relative contributions. Because the convective plumes have both ascending and descending 
motions, the average vertical mass fiux filters out their contribution and is therefore a good proxy of the 
amplitude of the acoustic modes. The left panel of Fig. [2] therefore displays the temporal evolution of pUz 
for the two simulations discussed before. An oscillatory behaviour is observed in both cases due to the radial 
oscillations of the fundamental acoustic mode excited by /^-mechanism. In the G8 simulation, the amplitude first 
grows exponentially until reaching the nonlinear saturation regime. At first glance, this time evolution looks very 
similar to what has been already observed in our former purely radiative simulations done in Gastine fc DintransI 
([2008b), that is, a linear growth of the amplitude and a saturation at a well-defined value. In contrast, the 
dynamics of the G8H8 simulation differ radically from the previous one as the amplitude remains weak and is 
highly modulated over time. No clear nonlinear saturation is observed in this case, meaning that the acoustic 
oscillations are more infiuenced by convective motions than in the previous DNS. 

To separate more precisely the relative contributions of the acoustic modes and the convective motions to the 
energy budget, t he velocity field of e ach simulations is projected onto an acoustic subspace built from normal 
eigenmodes (see Bogdan et al. 19931 ). Thanks to this formalism, it is possible to extract the time evolution 



of the kinetic energy contained in each acoustic mode found in the nonlinear simulation. The right panel of 
Fig. [2] displays the time evolution of the energy contained in the fundamental acoustic mode normalised by the 
total kinetic energy. In the G8 simulation, the acoustic energy linearly increases until its nonlinear saturation, 
in a similar way to what has been previously observed with the time evolution of pUz- Once this saturation 
is reached, 70% of the kinetic energy is contained in the radial oscillations of the fundamental acoustic mode, 
while the remaining is in the convective plumes. In other words, the acoustic oscillations are not affected much 
by the convective motions in this simulation. In contrast, this acoustic energy ratio remains very weak in the 
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Fig. 2. Left panel: Temporal evolution of the mean vertical mass flux pUz for the two simulations G8 (solid blue line) 
and G8H8 (solid green line). The two vertical dashed black lines define the boundaries of the zoom displayed in the 
bottom left corner. Right panel: temporal evolution of the energy contained in acoustic modes normalised by the total 
kinetic energy for the two simulations G8 (solid blue line) and G8H8 (dashed green line). 



other simulation. Despite some transient increases during which non-trifling values {c^ 10%) are obtained, the 
average ratio is less than 5%, and convective motions contains the bulk of the kinetic energy. In this case, 
the radial oscillations excited by the mechanism are thus quenched by convective plumes. This situation is 
relevant to the physics of Cepheids close to the red edge of the instability strip, where the unstable acoustic 
modes are supposed to be damped by the surface convective motions. This convective quenching of the acoustic 
oscillations may be the direct signature of the different density contrasts in the 08 and the G8H8 simulations: 
in fact, weaker stratification (as in G8H8) leads to bigger vortices (see Fig. [TJd), meaning that the energy is 
contained in larger convective structures. In our DNS, the am plitude of the /^-mechanism seems to be controlled 
by the screening effect due to these large convective vortices (^Gastine fc DintransI 



3 Limits of time-dependent convection models 



The nonlinear simulations of the convection-pulsation coupling, where the acoustic modes strongly modulate 
the convective motions over time (as in the G8 simulation) are also good candidates to test and compare the 
relevance of different prescriptions of 1-D time-dependent convection (TDC) models. We focus h ere on two 
popular formu lations widely used in Cepheids models, namely the TDC models of IStellingwer j ( 1982 ) and 
Kuhfui^l (|l986l ). In these formulations, a single equation for the turbulent kinetic energy St is added to the 
classical mean-field equations and the main second-order correlations, such as the convective flux, are expressed 
as a function of Et only: 



A 



Tstiz.t) = ast-£tsign{\/ - Vad)\/|V - Vad| 

ID 



(3.1) 



J^Kuiz, t) = aKuA\/ £t (V - Vad) , 

where V = dlnT/dlnp, Vad = 1 — Cv/cp and 



, A = cp {p) (T) and B = Jcp (T) Vad 



(3.2) 



where p is the pressure, p the density, and the brackets denote a horizontal average. These two TDC expressions 
i nvolves one dimensionless parameter alpha {ast and axu? respectively) that is poorly constrained by theory 
( Yecko et al.|[l998 ). The nonlinear results of the G8 simulation allow us to constraint it by using the x^-statistics 
so that to obtain its optimum value for each formulation. 

Figure [3] compares the convective f lux in the DNS with the tw o TDC prescriptions computed with the best- 
fit values obtained for ast and a-^u (jGastine fc DintransI 1201 Ibh . Stellingwerf's formulation seems to give a 
better agreement with the nonlinear simulation than KuhfuB's one. Indeed, the KuhfuB model overestimates 
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Fig. 3. Mean convective flux in the G8 simulation (solid black line), compared with the best TDC predictions based on 
the models of Stellingwerf (dashed blue line) and KuhfuB (dotted green line). 

the overshooting as the (negative) convective flux remains non-negligible until the bottom of the radiative zone. 
In contrast, the Stellingwerf proflle accounts for the local penetration of convective plumes better and shows the 
expected exponential-like decay in the negative convective flux when sinking in the radiative zone. However, the 
two models are fairly similar in the bulk of the convective zone, where convection is fully developed. One also 
notes that they both predict a negative flux at the top of the convective zone, which is an upper overshooting 
of convective motions near the surface that is not observed in the DNS. 

4 Conclusion 

The main weakness of all theories of 1-D time-dependent convection lies in the large number of free parameters 
involved in the description of convection. New constraints must therefore be found to reduce the intrinsic degen- 
eracy of these models and to check the relevance of the different assumptions underlying these parametrisations. 

Nonlinear 2-D direct numerical simulations are a useful way to address the convection-pulsation coupling 
that occurs in cold Cepheids. These simulations show a variable influence of convection onto the acoustic modes 
excited by mechanism: {i) either the amplitude of the acoustic modes remains very weak and convective 
motions quench the oscillations ; (n) or the kinetic energy is mainly contained in the acoustic modes and 
convective plumes are strongly modulated over time by the radial oscillations. While the former situation is 
relevant to the stabilisation of the oscillations of Cepheids close to their red edge, the latter is a good candidate 
to draw the limits of current TDC recipes. Focusing on two such widely used models, Stellingwerf 's formulation 
is found to give a better agreement with the nonlinear results than does KuhfuB's. 

This flrst comparison of TDC theories with fully nonlinear results coming from 2-D direct simulations of 
the convection-pulsation coupling emphasises how DNS can be helpful to validate and improve 1-D models of 
time-dependent convection. 

This work was granted access to the HPC resources of CALMIP under the allocation 2010-P1021 ( |http: //www. calmip. cict . fr[ ). 
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